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1 Introduction 



The workshop "HadAtom02" took place at CERN on October 14-15, 2002. It was the 4th in a 
series of workshops on bound states, in particular hadronic atoms. The previous ones were held 
in Dubna, Russia (May 1998) [1] and in Bern (October 1999, October 2001) [2,3]. The meeting 
was attended by about 40 physicists, and contributions were presented by 17 participants. 
The workshop programme included: 

• Hadronic atoms, in particular their • Experiments 

Production DIRAC at CERN 

Interaction with matter DEAR at DAFNE 

Energy levels PSI (Pionic Hydrogen Collaboration) 

Decays NA48 

• Meson-meson and meson-baryon scattering Others 

• K i4: decays 

The talks were devoted to recent experimental and theoretical progresses in the investi- 
gations of hadronic atoms. Among the highlights of the workshop were the presentation of 
preliminary results from the DIRAC collaboration on the measurement of the lifetime of pio- 
nium, as well as the first measurement of kaonic nitrogen by the DEAR collaboration. 

The speakers have provided a two-page summary of the presentations including a list of the 
most relevant references — these contributions are collected below. These proceedings include 
also a list of the participants with their e-mail, and the scientific program of the workshop. 
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Need for new low-energy tttt scattering data 



L. Girlanda 

ECT and INFN, Strada delle Tabarelle 286, 38050 Trento, Italy 

The importance of low-energy inr scattering data for probing the mechanism of spontaneous 
breaking of chiral symmetry (SBChS) in QCD has often been emphasized. Explicit calculations 
of the 7T7T scattering amplitude at two-loop level in the framework of chiral perturbation theory 
(ChPT) allow to relate the scattering parameters to the low-energy constants, describing the 
features of the vacuum. Present knowledge about these constants, complemented with standard 
assumptions about the size of the quark condensate (i.e. the validity of the Gell-Mann-Oakes- 
Renner relation) and available experimental information at medium and high energy, leads to 
very precise predictions for the ,S-wave scattering lengths [1] 

a° = 0.220 ±0.005, a 2 = -0.0444 ± 0.0010. (1) 

On the other hand, thanks to the solutions of the Roy Equations provided by Ananthanarayan, 
Colangelo, Gasser and Leutwyler [2], the scattering amplitude is parametrized in terms of only 
two parameters, the two S-wave scattering lengths, within very small uncertainty in the whole 
energy range below 800 MeV. Therefore any experimental information below 800 MeV can in 
principle be translated into a determination of the two S-wave scattering lengths, thus providing 
a test of the predictions of ChPT. We applied this procedure [3] to the recently published E865 
data [4] on charged K e 4 decays, supplementing them with available data in the 1 = 2 channel 
below 800 MeV. 

The ellipses in the figure (full lines) represent the 1-a and 2-a contours of our fit. Our 
analysis is exclusively based on the new solutions of the Roy Equations [2] and on direct ex- 
perimental information. In particular no use is made of the correlation between a® and a 2 , as 
inferred from the scalar radius of the pion [5]. If the dispersive determination of the latter, 
(r 2 ) s = (0.61 ± 0.04) fm 2 is used, this correlation, at the level of two-loop accuracy, is shown 
by the narrow curved band of the figure. The use of such a strong constraint would result 
in the filled shaded ellipse, which is perfectly compatible with the two-loop ChPT prediction, 
represented by the small filled ellipse. The two procedures yield different results at the 1-a 
significancy. However the narrow strip constraint depends substantially on the NNLO contri- 
butions, which involve unknown 0(p & ) low-energy constants: the same correlation without the 
inclusion of the 0(p 6 ) order would correspond to the narrow straight band in the figure. In 
order to obtain the 0(p 6 ) narrow strip one has to rely on estimates of the constants appearing 
at this order, based on the resonance saturation assumption. By reversing the argument, one 
can use the model-independent result for the scalar scattering lengths 

a° = 0.228 ±0.012, a 2 = -0.0382 ± 0.0038, (2) 

to test the narrow strip prediction, i.e. to determine the 0(p 6 ) constants. The latter are found 
unexpectedly large, which can be interpreted as a manifestation of the exceptional status of the 
scalar channel, characterized by a strong tttt continuum and OZI rule violation. The situation 
will hopefully be further clarified by the awaited results of DIRAC and the forthcoming K e4: 
experiment NA48-II. Unfortunately it is unlikely that the measurement of the pionium lifetime 
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can distinguish between the two procedures, since all the ellipses shown in the figure correspond 
to very similar values for the combination a[j — ctg. 

The fit results for a(] and can be used to extract the values of the subthreshold param- 
eters a and (3, of the low-energy constants Z3 and Z4 as well as of the main two-flavour order 
parameters: (uu) and F v in the limit m u = rrid = taken at the physical value of the strange 
quark mass. The results for the order parameters read 

X(2) = ^00 = 0.81 ±0.09, Z(2)=~ = 0.90 ±0.03, (3) 

signalling that the two-flavour standard ChPT is the appropriate expansion scheme. 

[1] G. Colangelo, J. Gasser and H. Leutwyler, Phys. Lett. B 488 (2000) 261 [arXiv:hep- 
ph/0007112]. 

[2] B. Ananthanarayan, G. Colangelo, J. Gasser and H. Leutwyler, Phys. Rept. 353 (2001) 207 
[hep-ph/0005297]. 

[3] S. Descotes, N. H. Fuchs, L. Girlanda and J. Stern, Eur. Phys. J. C24 (2002) 469 [hep- 
ph/0112088]. 

[4] S. Pislak et al. [BNL-E865 Collaboration], Phys. Rev. Lett. 87 (2001) 221801 [hep- 
ex/0106071]. 

[5] G. Colangelo, J. Gasser and H. Leutwyler, Nucl. Phys. B603 (2001) 125 [hep-ph/0103088]; 
Phys. Rev. Lett. 86 (2001) 5008 [hep-ph/0103063]. 
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Chiral perturbation theory and massive strange quark 

pairs 

Jan Stern 

Groupe de Physique Theorique, Institut de Physique Nucleaire, Universite de Paris- Sud , 

91 406 ORSAY - CEDEX 

Low - energy -kit scattering merely involves valence u and d quarks and it probes QCD 
in the vicinity of the SU(2) x SU(2) chiral limit in which m u = rrid = and m s is fixed 
at its physical value. Nevertheless, the influence of massive virtual ss pairs on low-energy 7r7r 
scattering, on the two-flavour chiral order parameters E(2) = — \im mutmd ^ (uu) \ ms =phys — F 2 B, 
F 2 = \im mu rnd ^ F 2 , as well as their impact on the estimates of low energy constants l 3 and U [1] 
may be important for the following two reasons: i) First,m s <C Ah ~ ^F w ~ lGeV suggesting 
some overall convergence of the SU(3) x SU (3) chiral expansion . On the other hand, as long as 
m s ~ Aqcd, the scale at which QCD becomes strong, the vacuum gets populated by abundant 
virtual ss pairs, whereas the creation of heavy quark pairs is suppressed, ii) The second reason 
is the existence of a strong long range correlation between strange and non-strange qq pairs in 
the J p = + channel: 

Z s = lim i [ dx(Tss(x)uu(0)) c (4) 

is subleading in large N c expansion and it is expected to be suppressed by the OZI - rule. 
Instead, recent sum rule estimates [2] suggest that Z s is of a normal size Z s ~ Aq CD , reflecting 
the violation of the OZI- rule and of the 1/N C expansion observed in the scalar channel. This 
phenomenological information was not available in the past and most of the original 
considerations and estimates [1], (including the G%PT variant) were based on the OZI-rule 
and large N c wisdom. Today, the strong vacuum OZI-rule violating correlation Z s provides a 
new element in the discussion and it naturally fits to our microscopic understanding of %SB 
in Euclidean QCD: The quark condensate represents the average and the correlator Z s the 
fluctuations of the density of smallest modes of the Dirac operator. Notice that Z s > 0. 

Vacuum fluctuations of ss pairs may lead to important m s dependence of the two-flavour 
condensate S(2). At m s = 0, S(2) coincides with the three-flavour condensate S(3) = 
— (uu) \ mu=md=m3=0 = FqB and its derivatives with respect to m s is precisely the correlator Z s 
(related to the low - energy constant L 6 ). Up to higher order terms one has 

S(2) = S(3)+m s Z s + .... (5) 

S(3) is the "genuine condensate" of QCD with three massless quarks and no massive quarks left 
which would be enough light to pollute the vacuum. m s Z s represents the positive contribution 
to the two-flavour condensate S(2) that is induced from massive ss vacuum pairs. The two 
contributions can be of sensibly comparable size [2], reflecting the suppression of E(3) relative 
to S(2) due to vacuum fluctuations. This can be seen from a non-perturbative analysis of 
Goldstone-boson mass and decay constant Ward identities [3], provided the quark mass ratio 
r = m s /m is not too small (r > 15). For small Z s , which in turn implies a precise fine tuning 
of L 6 (M P ) to the critical value L^ rit = -0.26 x 1(T 3 , one would have S(2) w E(3). For L 6 
slightly above Lg"*, the fluctuations grow and the yield of the genuine condensate S(3) in 
Eq. ©decreases. This decrease is, however, compensated by a growing induced contribution 
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m s Z s . As a result, the massive ss pairs stabilize the two-flavour condensate S(2) and keep it 
large independently of vacuum fluctuations and of the suppression of the genuine condensate 
E(3) ,(see the talk by S. Descotes-Genon). 

It is instructive to compare the Nf = 2 %PT (i.e. expansion in powers of m u = = m for 
a fixed m s ) with the Nf = 3 %PT expansion (in powers of m and m s ) of the same quantity 
F%M%. Keeping the LO and the NNLO contributions, one gets 

FlMl = 2mS(2) + ^(4Z 4 - l 3 ) + . . . (6) 

FlMl = 2mS(3) + 2mm s Z s + Am 2 (Z s + A) + . . . . (7) 

respectively. The two expansions have to coincide order by order in m. The most important 
part of the OZI - rule violating contribution to rhs of Eq. (JZJ), m s Z s , is absorbed into S(2). 
The contribution of Z s to the LEC's Z3, Z4 is not enhanced by the factor m s and it should be of 
a comparable size as the "normal" contribution A (related to L 8 ). Hence, vacuum fluctuations 
will slightly shift the OZI-rule based estimate [1] Z3 = 2.9 ± 2.4 towards a negative values 
without spoiling the convergence of the Nf = 2 x?T. The Nf = 2 expansion is dominated by 
the condensate term E(2), since the latter is enhanced by the contribution induced by ss pairs. 
On the other hand, the genuine condensate S(3) need not dominate the Nf = 3 expansion of 
Goldstone boson masses, since the vacuum fluctuation term m s Z s now counts as NLO. It is 
useful to express both condensates in appropriate GOR units 

-WH 2 ^ ", = 2,3 (8) 

7T 7T 

and to compare X(2) > X(3) with one. The Nf = 2 low-energy parameters can be extracted 
from precise irn scattering experiments . Our analysis [4] including most recent experimental 
data yields X(2) = 0.81 ± 0.07 , F 2 /F* = 0.89 ± 0.02 and l 3 = -17 ± 15 confirming both 
the dominance of the two-flavour condensate and the expected shift of /3. (See the talk by L. 
Girlanda for more details.) 

Suppression of X(3) does not necessarily mean a bad convergence of the Nf = 3 for 
"good observables" such as ((7|). However, provided in the latter the LO term 2mE(3) and the 
NLO term 2mm s Z s are of a comparable size, we face an instability which prevents us to proceed 
perturbatively eliminating low-energy parameters in terms of observables. A systematic non- 
perturbative alternative to this standard procedure does exist and it makes a new Nf = 3 
analysis of precise 7T7T , 7rK and ?7-decay data both meaningful and interesting. 

[1] J. Gasser and H. Leutwyler, Annals Phys.158 (1984) 142, Nucl.Phys. B250 (1985), 465. 
[2] B. Moussallam , Eur. Phys. J. C14 (2000) 111, [ arXiv:hep-ph/9909292], S. Descotes, JHEP 

0103 (2001) 002, [ arXiv: hep-ph/0012221]. 
[3] S. Descotes-Genon, L. Girlanda, and J. Stern, [arXiv: hep-ph/0207337] and references 

therein. 

[4] S. Descotes, N.H. Fuchs, L. Girlanda and J. Stern, Eur.Phys.J. C24 (2002) 469, [arXiv: 
hep-ph/0112088]. 
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QED Theory of Hydrogen-like Atoms 

Savely G. Karshenboim 

D. I. Mendeleev Institute for Metrology (VNIIM), St. Petersburg 198005, Russia 
Max-Planck-Institut fur Quantenoptik, 85748 Garching, Germany 

About fifty years ago the experimental discovery of the Lamb shift and the anomalous 
magnetic moment of electron initiated a long competition between QED theory and experiment. 
Today a list of quantities used to test QED is quite impressive and includes: 

• the anomalous magnetic moments of the electron and muon; 

• hyperfine splitting of the Is and 2s states in hydrogen and deuterium atoms and helium-3 
ion; 

• gross and fine structure intervals and the Lamb shift in hydrogen, deuterium and helium-4 
ion; 

• hyperfine splitting and Is — 2s transition in muonium and positronium, 

• decay rates of para- and orthopositronium including rare modes; 

• energy levels of high- (like e.g. bismuth) and medium-Z (like e.g. carbon) ions with few 
electrons; 

• g factors of a bound electron in hydrogen-like atoms; 

• energy levels of muonic and exotic atoms; 

• energy levels of three-body atoms (helium, antiprotonic helium, muonic helium). 

Twenty years ago the experimental accuracy of most QED tests was significantly better 
than the theoretical one. With time situation changed so that after successful devopments 
modern QED theory dominates over experiment. 

However, there is the third side in this competition. Unfortunately, theory and experiment 
cannot speak the same language. Experiment leads to results expressed in units of Hertz and 
Mev, while theory delivers only expressions containing values of fundamental constants (such 
as the fine structure constant a, electron mass m e ) and particle/nuclear parameters (like e.g. 
the proton charge radius). In other words, theory is actually not in position to predict any 
numbers building instead bridges between different experiments. As a bridge, theory is the 
most accurate part of most QED tests. However, trying to calculate some quantity needed to 
make a prediction we often discover that the input data of the calculation, i.e. the fundamental 
and auxiliary constants, are not known accurately enough from their measurements in other 
experiments. 

In the case of hadronic effects for some QED quantities (like e.g. hyperfine interval in the 
hydrogen atom or anomalous magnetic moment of muon) we even need functions (such as the 
proton electric form factor) to be somehow determined as the input data. 

The lesson we have to learn from precision studies of simple atoms and QED tests is that the 
situation for free particles (g — 2 for electron and muon), low Z atoms (hydrogen, muonium), 
high Z atom (bismuth, lead), muonic and exotic atoms is actually quite similar: 
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• the QED calculations themselves can be done with high enough accuracy to compete with 
experiment. Actually, QED theory rather wins the competition. 

• However, that is not enough: 

— we need to have access to accurate values of the fundamental constants to transform 
QED expressions into numbers; 

— we also need to know some functions like form factors to put into the calculation. 

• Eventually, comparison of theory and experiment is often not quite a QED test, but rather 
a determination of some fundamental constants or study of nuclear effects. 

The crucial orders of magnitude important for comparison of theory and experiment are 
summarized in the Table below. In most of the cases the QED tests confirm theory. 



Quantity Order 

Hydrogen, deuterium (gross structure) a(Za) 7 m, a 2 (Za) 6 m 

Hydrogen (fine structure and Lamb shift) a(Za) 7 m, a 2 (Za) 6 m 

3 He + ion (Is /2s hyperfine structure) a(Za) 7 m 2 /M,a(Za) 6 m 3 /M 2 , 

a 2 (Zafm 2 /M, (Za) 7 m 3 /M 2 

4 He + ion (gross structure and Lamb shift) a(Za) 7 m, a 2 (Za) 6 m 

Muonium (Is hyperfine structure) (Za) 7 m 3 /M 2 , a(Za) 6 m 3 jM 2 , 

a(Za) 7 m 2 /M 

Positronium (gross structure and Is hyperfine structure) a 7 m 

Parapositronium (decay rate) a 7 m 

Orthopositronium (decay rate) a 8 m 

Parapositronium (47 branching) a 8 m 

Orthopositronium (57 branching) a 8 m 

Table 1. Comparison of QED theory and experiment: crucial orders of magnitude 

More detail on comparison of theory versus experiment can be found in [1]. A detailed 
review of theory is presented in [2]. Most of recent results and projects related to QED tests 
and simple atoms along with reviews were presented at two conferences on Precision Physics of 
Simple Atomic Systems (PSAS) in 2000 (Castiglione della Pescaia) [3] and 2002 (St. Petersburg) 
[4]- 

This work was supported in part by the RFBR grant 00-02-16718. 

[1] S. G. Karshenboim, In Atomic Physics 17 (AIP conference proceedings 551) Ed. by E. 

Arimondo et al. (AIP, 2001), p. 238; arXiv:hep-ph/0007278. 
[2] M. I. Eides, H. Grotch and V. A Shelyuto, Phys. Rep. 342, 63 (2001). 
[3] S. G. Karshenboim, F. S. Pavone, F. Bassani, M. Inguscio and T. W. Hansch. Hydrogen 

atom: Precision physics of simple atomic systems (Springer, Berlin, Heidelberg, 2001). 
[4] The contributed papers to PSAS 2002 are scheduled for publication in Canadian Journal 

of Physics in November, 2002, while the book of review contributions is in preparation to 

be published in 2003. 
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Ground-state energy of pionic hydrogen to one loop 



J. Gasser 

Institute for Theoretical Physics, University of Bern, Sidlerstrasse 5, CH-3012 Bern 

In my talk I discussed recent work on pionic hydrogen performed in Ref. [1]. In Ref. [2], 
the relation between the scattering length combination \a^ + — a^ + \ and the strong energy-level 
shift of pionic ground state has been worked out at leading order in the low-energy expansion 
of the isospin breaking correction. In Ref. [1], we have carried out the calculation at next- 
to-leading order, where the effect of the photon-neutron intermediate state first shows up in 
the chiral expansion. This amounts to evaluating the elastic n~p — > ir~p amplitude at order 
p 3 , including isospin breaking contributions generated by the quark mass difference m u — 
and by virtual photon loops. In order to have a systematic framework, we have extended the 
infrared regularization method developed by Becher and Leutwyler [3], such that it includes 
virtual photons as well. Furthermore, we have performed a heat-kernel evaluation of all singular 
parts in the one-loop calculation, including meson, baryon and photon fields. As far as we are 
aware, this provides for the first time a consistent set of effective lagrangians at this order in the 
low-energy expansion. We find that triangle-type graphs generate large corrections, whereas 
the effect of the low-energy constants at order p 3 is suppressed by the factor rrii T /m pTOton . On 
the other hand, the constant /i, that occurs at leading order, may contribute significantly to 
the energy shift [1]. Unfortunately, it has not yet been determined with sufficient accuracy. We 
conclude that the potential-model calculations available in the literature largely underestimate 
the systematic error, because these models are not properly matched to the underlying theory 
(QCD+QED). Furthermore, claims made in the literature that a precise measurement of the 
strong energy-level shift in pionic hydrogen will provide a corresponding precise determination 
of I too optimistic at the present stage of the theoretical knowledge: Experiment 

is ahead of theory in this case. For further references and details, I refer to [1]. A very recent 
article is [4] - we plan to comment on this work elsewhere [5]. 

[1] J. Gasser, M. A. Ivanov, E. Lipartia, M. Mojzis and A. Rusetsky, Eur. Phys. J. C 26 (2002) 

13 [arXiv:hep-ph/0206068]. 
[2] V.E. Lyubovitskij and A. Rusetsky, Phys.Lett. B494 (2000) 9 [arXiv:hep-ph/0009206]. 
[3] T. Becher and H. Leutwyler, Eur. Phys. J. C9 (1999) 643 [arXiv:hep-ph/9901384]. 
[4] T.E.O. Ericson, B. Loiseau and S. Wycech, arXiv:hep-ph/0211433. 
[5] J. Gasser and A. Rusetsky, in preparation. 
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The influence of the channel on the energy and width 

of pionic hydrogen 

G.C. Oades 1 , G. Rasche 2 and W.S. Woolcock 3 

1 Institute of Physics and Astronomy, Aarhus University, DK-8000 Aarhus C, Denmark. 
2 Institute fur Theoretische Physik der Universitat, Winterthurerstrasse 190, CH-8057 Zurich, 

Switzerland. 

3 Department of Theoretical Physics, IAS, The Australian National University, Canberra, 

ACT 0200, Australia. 

We reconsider the corrections presented in the paper of Sigg et al. [1] for the extraction 
of the hadronic s— wave 7rN scattering lengths from measurements of the energy and width 
of pionic hydrogen. Our starting point is the Deser type formulae given in the multi-channel 
situation by Rasche and Woolcock [2]. These relate the energy shift to the nuclear scattering 
length, a" c , for the process ir~p — > ir~p and the width for decay to the 7T°n channel to the 
scattering length, a™ for the process ir~p — > n°n. These nuclear scattering lengths are in turn 
related to the hadronic scattering lengths by the expressions 

c& = c&(l + 

a n c0 = a h c0 (l + 5 c0 ), 

where 5 CC and 5 c q are the electromagnetic corrections. Values for these were first considered in 
[2], later examined in more detail in [1] and more recently recalculated in [3]. 

In addition to the hadronic channels n~p and 7r°n, pionic hydrogen also decays to the 
channel 771, the widths being related by the Panofsky ratio 

^ = P = 1.546 ±0.009. 

l-\n 

This channel also influences the EM corrections and the first crude estimates of its effect were 
made in [1] where the effective potential in the 7r~p — > ir°n channel was empirically increased 
so that it gave a value a e JJ which reproduced the full width i.e. 

This modified potential was then used to calculate new EM corrections and the differences 
between these values and the values from the unmodified 2-channel calculation, A<5 CC and A5 c0 , 
were assigned to the influence of the •yn channel. The final values with errors given in [1] were 

5 CC = -(2.1 ±0.5)% 
5 c0 = -(1.3 ±0.5)% 

where the contributions from the yn channel were 

A5 CC = -0.73% 
A5 c0 = +0.15%. 
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What the authors of [1] overlooked was that the modification of the effective potential in the 
7i~p — > 7r°n channel not only changes a™ to aJJ but also modifies a" c , an effect which can also 
be considered as a type of EM correction. In order to examine this effect we have extended the 
calculations used in [3] to the case where we have the three channels, 7r~p, Tr°n and •yn. Such 
calculations involving mass photons are possibly dangerous due to the fact that the reduced 
mass of the yn system is 0. As pointed out in [4], this problem does not occur in the case of a 
relativistic equation of motion. In our case we use the relativised Schrodinger equation where 
the appropriate term in the non-relativistic equation, 2M^ d V CJ , becomes 2M^ d f cl V n where 

ra Jn M n + m 7 2m y W 

which has a finite limit as m 7 0. In this calculation we also need effective potentials to 
describe the three channel system. For the n~p and 7r°n channels we use hadronic potentials 
fitted to the low energy phase shifts together with the Coulomb and vacuum polarization 
potentials corrected for finite size effects. For n~p — > yn and 7r°n — > yn we adjust the potentials 
to fit the E 0+ photoproduction multipole amplitudes and for yn — > yn we adjust to fit the low 
energy yn — > yn cross section. We now make the same type of calculations as in the two 
channel case i.e. we compare the value of the threshold K- matrix calculated in the full three 
channel case with its purely hadronic value. For the EM corrections, preliminary calculations 
show only a very small change in 5 c q but 5 CC increases to a value around zero. 

Can we understand this physically? What we are seeing might well be called an extended 
Ball-Frazer mechanism. The addition of the extra inelastic channel gives an effective attraction 
which increases the binding energy. This seems to be an effect which comes from the require- 
ments of multichannel unitarity and as such will not be seen in lowest order calculations which 
ignore these requirements. 

Another very important point is that our EM corrections refer to a hypothetical hadronic 
world which is, by tradition, defined as one in which nucleons have the mass of the proton, 
pions the mass of the charged pion and where the effective potentials are isopin invariant. 
This should be remembered when comparing to the results of Gasser et al. [5] who refer to a 
hadronic situation consistent with chiral perturbation theory. We could, of course, use a similar 
definition of our hadronic reference point. The hadronic masses would then be taken from chiral 
perturbation theory and the effective potentials adjusted to reproduce the low energy chiral 
perturbation theory amplitudes. To do this without assuming isospin invariance requires the 
low energy amplitudes in all three channels, n~p — > n~p, n~p — > n°n and 7i°n — > 7i°n and we 
therefore ask the chiral perturbation theory experts for these values so that we can proceed 
with this work. 

[1] D. Sigg, A. Badertscher, P.F.A. Goudsmit, H.J. Leisi and G.C. Oades, Nucl. Phys. A 609 
(1996) 310. 

[2] G. Rasche and W.S. Woolcock, Nucl. Phys. A 381 (1982) 405. 

[3] A. Gashi, E. Matsinos, G. Rasche, G.C. Oades and W.S. Woolcock, arXiv:hep-ph/0009081 

and Nucl.Phys A 686 (2001) 463. 
[4] PB. Siegel and W.R. Gibbs, Phys. Rev. C 33 (1986) 1407 and W.B. Kaufmann and 

W.R. Gibbs, Phys. Rev. C 35 (1987) 838. 
[5] J. Gasser, M.A. Ivanov, E. Lipartia, M. Mojzis and A. Rusetsky, arXiv:hep-ph/0206068. 
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Coupled channel approach to breakup of pionium 



T.A. Heim 1 , K. Hencken 1 , M. Schumann 1 , D. Trautmann 1 and G. Baur 2 

1 Institut fur Physik, Universitat Basel, Klingelbergstrasse 82, CH-4056, Basel, Switzerland 
2 Institut fur Kernphysik, Forschungszentrum Julich, 52425 Jiilich, Germany 

Since the experiment DIRAC does not directly observe the neutral pions from the strong decay 
channel 7r + 7r~ — > 2tt°, but rather looks at the charged pions from electromagnetic breakup in 
the target, this electromagnetic breakup process must be calculated with very high accuracy 
(on the order of 1%) as a required input for the analysis of the experiment. In order to extract 
the pionium lifetime from the observed breakup probability, the passage of pionium through 
the target matter is studied in a simulation which in turn requires detailed knowledge of the 
cross sections for bound-bound and bound-free transitions of the pionium [1]. 

In previous work [2-4] we have shown that the semiclassical approximation provides the 
appropriate tool for this problem. In first order Born approximation the accurate description of 
both the pionium and the target system with its atomic structure, as well as their relativistically 
correct interaction (including magnetic terms), has been implemented successfully (accurate to 
far better than 1% within this model). However, it has been noted by studying the total cross 
sections [5], and veryfied by explicitly calculating the bound-bound transitions in the Glauber 
formalism [6] that higher order contributions can change the first order results by up to several 
per cent. We have estimated the validity of the Glauber approximation in [6] by studying the 
influence of finite interaction times within a simple extension of the sudden approximation. We 
found that the results varied within less than 0.1% which may thus be regarded as the accuracy 
of the higher order calculation (clearly below the required 1% limit). 

In order to provide a further, independent test 
of the validity of the Glauber approximation, 
we study the time evolution of the occupation 
probabilities for a few selected bound states 
of the pionium in a coupled channel approach. 1 
The figure shows this time evolution for various e 
channels as indicated. In this example, the pio- J 
nium starts out initially in the 3s-state (all other In- 
states not occupied; Nickel target). As time ° 
evolves, other pionium states acquire non-zero 
occupation probabilities. Our goal is to verify 
whether the cross sections obtained in this ap- 
proach agree with our previous (numerically less 
demanding) calculations in a few test cases. 




v t [atomic units] 



[1] L.G. Afanasyev et al, Phys. Lett. B 338 (1994) 478. 

[2] Z. Halabuka et al, Nucl. Phys. B 554 (1999) 86. 

[3] T.A. Heim, K. Hencken, D. Trautmann and G. Baur, J. Phys. B 33 (2000) 3583. 

[4] T.A. Heim, K. Hencken, D. Trautmann and G. Baur, J. Phys. B 34 (2001) 3763. 

[5] L. Afanasyev, A. Tarasov and O. Voskresenskaya, J. Phys. G 25 (1999) B7. 

[6] M. Schumann et al, J. Phys. B 35 (2002) 2683. 
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On the influence of low n and high n states cross section 
in the break-up probability of pionium 

C. Santamarina 1 and L.G. Afanasyev 2 

1 Institute of Physics, University of Basel, CH-^056, Basel Switzerland 
2 Joint Institute for Nuclear Research, Dubna, Moscow Region, 141980 Russia 

The study of the break-up probability (P&r) dependence of pionium as a function of lifetime 
in a target of definite material and thickness is a master key in the measurement that DIRAC 
collaboration aims to achieve in its experiment PS-212 at CERN [1]. 

A detailed explanation of the physical problem can be found in [2] and [3] and can be resumed 
as to solve the differential equation system that links the population p n i m of the different atomic 
bound states (defined by the usual quantum numbers of an hydrogen-like system, n, I and m) 
as a function of the target position of the atom (s) : 



dpnlm(s) s—*. n'l'm' ( \ 

= 1^ a nlm Pn'l'm'{S) 



ds 



n'l'm' 



n'l'm' 
nlm 

-.nlm _ 
''nlm 



~total 



pN 



if / = 



2M^ 

PCT„ 

if I ^ 



where p is the density and A the atomic weight of the target, Nq the Avogadro number, M n 
the mass of the pion, r„ the lifetime of the bound state and P the center of mass momentum of 
the atom (average around 4.2 GeV/c in DIRAC experimental conditions). However, the main 
input to the equation is given by the pionium-target atom interaction cross sections, either 
for the probability of a transition between bound states (o"™/'™" 1 ') or for the total interaction 
probability «f^). 

These cross sections have been calculated in the Born approximation using different pa- 
rameterizations for the atomic form factors [2,4] and with the Glauber [5] approximation that 
accounts for multi-photon exchange. The three calculations lead to similar results with discrep- 
ancies between 1% and 4% if the states involved in the transition are low energy ones (n < 5) 
and increase as n does up to differences higher than 10%. 

However, the difference in the breakup probability result of pionium is to the level of 1%. 
As an example, asumming 2.91 femtoseconds for pionium lifetime: 



pG 



0.454 



pBl 



0.460 



pB2 



0.465 



where P£ was calculated with the Glauber cross sections of [5], P^ r 1 with the Born cross sections 
of [2] and P^ r 2 with the Bron cross sections of [4] . 

This agreement between the three results is explained by four things: 

1. The atomic bound states are produced (initial conditions) according to: p p n rod oc 1/n 3 [6] 
and hence the amount of atoms created with n > 5 is negligible. 



The probability of annihilation also behaves as p^ mh 
almost no atoms annihilated from states with n > 5. 



oc 1/n 3 [7] and hence there are 



3. The atoms show a clear tendency to be excited or broken rather than to be de-excited, 

J2n'<n,l',m' ^nllrT ^ J2n'>n,l' ,m' ^nllrT (see Figure). 
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4. The total cross section increases very fast as n does. And as a consequence the mean free 
path of the atom (X n i m = A/ a l °^ pN Q ) decreases very fast with n (see Figure). 

All this means that the states with n > 5 are populated by excitations from lower n states 
since very rarely atoms are directly created in these states. Moreover, de-excitations from these 
states to lower n states very seldom happen also. This means that the solution of the differential 
equation system for the low n states (where annihilation probability is not negligible) does not 
depend on the population of high n states and hence, the discrepancies between the different 
sets of cross sections will not lead to differences in the solutions of the equation system for these 
states populations. Meanwhile, an atom which is excited into a high n state has a very small 
free path (in all the three sets of cross sections) compared to the target dimensions of DIRAC 
(tenths of micrometers) and will interact with very high probability until it ends up broken. 
The discrepancies between the different sets of cross sections will be then irrelevant since they 
all lead to a breakup probability result of practically 100% for a bound state with n > 5 4 . 
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Figure 2: On the left we show up the comparison between the excitation and the de-excitation 
cross sections of [2] for the S sates of the atoms. On the right we can observe the mean free 
path value of S states. 

We want to thank Prof. Dr. Dirk Trautmann and his group for have kindly provided us 
their cross sections sets. 
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4 However this would not be the case if the target thickness was much smaller. 
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Recent results from DEAR at DA<E>NE 



C. Guaraldo on behalf of the DEAR Collaboration 
LNF-INFN, Via E. Fermi 40, OOO44 Frascati (Roma), Italy 

In May 2001 the DEAR (DA$NE Exotic Atom Research) collaboration [1] performed the 
first measurement of an exotic atom (kaonic nitrogen) at the DA$NE collider of the Laboratori 
Nazionali di Frascati dell'INFN. The 7 — > 6 at 4.6 keV and 6 — > 5 at 7.6 keV transitions were 
measured with 3 er and 3.7 a statistical significance, and the results published in [2]. After this 
first measurement, a second one was performed in May- June 2002, with an improved setup and 
better machine conditions (new optics which allowed the increase of luminosity and the decrease 
of background). The goals of this measurement were: to obtain a clean spectrum, in which 
only kaonic nitrogen lines are present, to reduce the continuous background with respect to the 
previous run (May 2001) in the view of the future kaonic hydrogen measurement and to perform 
a fine-adjustment of the degrader, which takes into account the <fi production mechanism at a 
machine as DA$NE (i.e. the boost). All these goals were successfuly reached. 

The target used for this measurement, differently from before, was done in kapton 75/im 
width, reinforced with carbon fiber. This had as positive consequence the fact that in the X 
ray spectrum all the eventual electronic transition lines corresponding to "biasing signals" , as 
iron, manganese, copper, were absent or much reduced, and the konic nitrogen signal was much 
clearer. The measurement was performed with the target filled with nitrogen at about 78 K and 
0.98 bar (p ~ 4.3p^Tp)- As detector, 16 CCD-55 were used. A pure-background spectrum (so- 
called "empty-target" ) was obtained as well, by using an over-dimensioned degrader, stopping 
the kaons inside, so not allowing the formation of the kaonic nitrogen atoms inside the nitrogen 
target. This spectrum was used for the shaping of the background and to obtain a background 
subtracted spectrum. 

A preliminary analysis of the data gave as result an overall X ray kaonic nitrogen spectrum 
which is shown in Figure 1, together with a fit in the region of interest. The empty-target 
subtracted spectrum is shown in Figure 2. The calcium electronic transition is due to the 
presence of calcium in the carbon fiber reinforcement of the target, while zirconium was placed, 
as a thin foil, for calibration purpose. 

The analysis (preliminary) of this spectrum gave a number of events in the 7.6 keV peak 
of 1400 ± 132 (10.6 a statistical significance) and in the 4.6 keV peak of 700±142 (about 5 a). 
More refined analysis are going to be performed. 

The spectrum in Figure 1 shows the nice feature of a continuous backgrouns; no electronic 
transitions are present in the region of interest (4-10 keV), apart eventually from some iron, 
which will be shielded in the next run. Moreover, a reduction of background by a factor 
about 5 was reached with respect to May 2001 run and the optimisation of the degrader was 
successfuly performed. The kaonic nitrogen transitions were identified with a much better 
statistical significance as well. 

The next stage of the experiment is the measurement of the K-complex in kaonic hydrogen 
transitions. 

[1] S. Bianco et al, Rivista del Nuovo Cimento 22, No. 11, (1999) 1 
[2] G. Beer et al, Phys. Lett. B535 (2002) 52. 
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Figure 3: The kaonic nitrogen X ray spectrum with a fit in the region of interest. 
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Figure 4: The kaonic nitrogen empty-target subtracted spectrum. 
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The Pionic-Hydrogen Experiment at PSI 



D. F. Anagnostopoulos 1 , S. Biri 2 , H. Fuhrmann 3 , D. Gotta 4 , M. Giersch 3 , 
A. Gruber 3 , A. Hirtl 3 , M. Hennebach 4 , P. Indelicato 5 , Y.-W. Liu 6 , B. Manil 5 , 
V. M. Markushin 6 , N. Nelms 7 , L. M. Simons 6 , P. A. Schmelzbach 6 , 
M. Trassinelli 5 , and J. Zmeskal 3 

1 Department of Material Science, University of Ioannina, GR-J^5H0 Ioannina, Greece, 

2 Institut of Nuclear Research; Hungarian Academy of Science, H-J^OOl Debrecen, Hungary 

3 IMEP, Osterreichische Akademie der Wissenschaften, A- 1090 Vienna, Austria, 4 Institut 
fur Kernphysik, Forschungszentrum Julich, D-52^25 Jiilich, 5 Laboratoire Kastler-Brossel, 
Universite Pierre et Marie Curie, F-75252 Paris, France, 6 Paul-Scherrer-Institut (PSI), 

CH-5232 Villigen, Switzerland, 7 Department of Physics and Astronomy, University of 

Leicester, Leicester LEI7RH, England 

In pionic hydrogen the hadronic pion-nucleon interaction manifests itself by a change of the 
energies and of the natural line width of X-ray lines as compared to a purely electromagnetically 
bound atomic system. Experimentally accessible are the transitions to the Is ground state 
emitted in the last de-excitation step of the atomic cascade. Any observed strong-interaction 
effect can be attributed fully to the Is state, because 2p-state effects are negligibly small. 
Hence, the s-level shift and width are exclusively owing to the pion-nucleon s-wave interaction, 
which is described in the limit of isospin symmetry by the isoscalar and isovector scattering 
lengths a + and a~ [1]. 

To improve on the accuracy for the hadronic parameters as compared to previous measure- 
ments [2], a thorough study of a possible influence of de-excitation processes is essential. A 
first series of measurements has been completed by the new pionic-hydrogen experiment at 
the Paul-Scherrer-Institut (PSI), using the new cyclotron trap, a cryogenic target and a Bragg 
crystal spectrometer equipped with spherically bent silicon and quartz crystals and a large-area 
CCD array [3]. Data analysis is in progress. 

In order to identify radiative de-excitation of the irH system - when bound into complex 
molecules formed during collisions n~p + H 2 — > [(ppn~)p}ee [4] - the energy of the nH(3p— Is) 
transition was measured at various target densities. X-ray transitions from molecular states 
should show up as low-energy satellites with density dependent intensities because of different 
collision probabilities. In our experiment, covering the pressure range from 3.5 bar to liquid, no 
density effect could be established [5] . It is concluded that the decay of molecules is dominated 
by Auger emission. The value derived for the hadronic shift is in agreement with the result of 
the previous experiment. 

At present, the accuracy for the hadronic broadening (7%) is limited by a not precisely 
known correction to the measured line width originating from the Doppler broadening due 
to Coulomb de-excitation [2]. For that reason the precisely measured ls-level shift in pionic 
deuterium was used together with the shift of hydrogen in the determination of the nN scat- 
tering length [6]. This procedure, however, requires a sophisticated treatment of the 3-body 
system ttD. In addition, up to now it cannot be excluded that the radiative decay channel after 
molecule formation is enhanced in deuterium compared to hydrogen. 

To improve the data base on the line broadening due to Coulomb de-excitation, the three 
nH(2p - Is) (2.4 keV), irH(3p - Is) (2.9 keV) and irH(4p - Is) (3.0 keV) transitions were 
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studied at a target density equivalent to 10 bar. An increase of the line width was found for the 
2p — Is line compared to the 3p — Is transition, which is attributed to the higher energy release 
available for the acceleration of the pionic-hydrogen system. This result is corroborated by a 
reduced line width of the Ap — Is line. For the total line width of the nH(3p — Is) transition an 
increase of 5-10% was found compared to the result of Schroder et al. [2], which may be due 
to the significantly improved background conditions in the new experiment. 

The response of the crystal spectrometer was obtained from the 7r 12 C(5g— 4/) line (3.0 keV), 
which is negligibly narrow compared to the experimental resolution (Fig. 1). The statistics of 
such a measurement, however, and with that the accuracy for the response function, is limited 
due to beam time considerations. In addition, the crystal response has been measured with 
X-rays emitted from helium-like Ar, ionised by means of an electron-cyclotron resonance ion 
trap (ECRIT) set up at PSI [7]. With that, studies of the Bragg crystals used up to now with 
good statistics could be performed within a reasonable time scale. 




energy/eV energy/eV 

Figure 5: Fig. 1: Ground-state transition 3p — Is in pionic hydrogen and the pionic carbon 
5 — 4 transitions measured with a quartz crystal. 

^From about 2005 on, Coulomb de-excitation will be studied in detail in the absence of 
strong-interaction effects by measuring K transitions from muonic hydrogen. Together with 
the detailed knowledge of the response function by using the ECRIT and a newly developed 
cascade code [8], which includes the velocity dependence of the atomic cascade, a sufficiently 
accurate correction for the Doppler broadening in pionic hydrogen should be achieved to extract 
the hadronic broadening at the level of about 1%. 

[1] V. E. Lyubovitzkij and A. Rusetski, Phys. Lett. B 494 (2000) 9. 
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Lifetime measurement of tt + 7t atom at DIRAC 

V. Yazkov on behalf of DIRAC Collaboration 

Skobeltsyn Institute for Nuclear Physics of Moscow State University 

Introduction 

Pionium or A 2 ^ is a hydrogen- like atom consisting of tt + and tt~~ mesons. The lifetime of 
this atom is inversely proportional to the squared difference between the S-wave irn scattering 
lengths for isospin and 2, |ao — 02 1- This value is predicted by chiral perturbation theory 
(ChPT) [1], and a measurement of the 7r + 7r~ atom lifetime provides a possibility to check 
predictions of ChPT in a model-independent way. 

Method of lifetime measurement 

The A 2n are produced by Coulomb interaction in the final state of -k + -k~ pairs generated in 
proton-target interactions [2,3]. After production A 2n travel through the target and some 
of them are broken up due to their interaction with matter: 'Atomic pairs" are produced, 
characterised by small pair cm. relative momenta Q < 3 MeV/c. These pairs are detected 
in the DIRAC setup. Other atoms annihilate into 7r°7r°. The amount of broken up atoms ua 
depends on the lifetime which defines the decay rate. Therefore, the breakup probability is 
a function of the A 2n lifetime. The dependence of Pbr on the lifetime r is determined by the 
solution of differential transport equations [4]. 

Also 7r + 7r~ pairs are generated in free state. Essential fraction of such pairs ("Coulomb 
pairs") are affected by Coulomb interaction, too. 

The aim of DIRAC is to measure the A 27T breakup probability P\x(t). -Pbr(T) is the ratio 
between the observed number of "atomic pairs" and the number of produced 7r + 7r~ atoms which 
is calculated from the measured number of "Coulomb pairs" . 

Experimental setup 

The purpose of the DIRAC setup [5] is to detect ir + 7i~ pairs with small relative momenta. This 
setup is located at the CERN T8 beam area (East Hall). It became operational at the end of 
1998 and uses the 24 GeV proton beam from PS accelerator. 

The setup resolution over the relative cm. momentum Q of 7r + 7r~ pair is better than 
1 MeV/c. 

Experimental data 

The experimental Q-distribution of tt + tt~ pairs is fit by an approximated distribution of "free 
pairs" in the region Q > 3.5 MeV/c where "atomic pairs" are absent. The number of "atomic 
pairs" is obtained as excess of experimental Q-distribution above the approximated distribution 
of "free pairs" only in the region Q < 3.5 MeV/c. Differences of experimental and approximated 
distributions are shown in Fig. Elfor data collected in 1999, 2000, 2001 and in 15 days of 2002. 
Platinum, nickel and titanium targets havebeen used. 
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Figure 6: Signals of "atomic pairs" in the region Q < 3 MeV/c are obtained with different 
targets. For 2002 only data collected in 15 days are analysed. 



The overall statistics of "atomic pairs" which is collected before the end of 2001 is more 
than 5000. It provides statistical accuracy for the lifetime measurement at the 20% level. 
The data analysis procedure is not completed yet and only very preliminary lifetime estima- 



tions are r = 2.8±H ■ 10 15 s for data collected in 2000 with nickel target and r = 5.4l ■ 10 15 s 



for all titanium target data. Averaged value is r = 3. 6^0!? • 10 15 s. 



Conclusions 

Preliminary results have been achieved by analysing a sample of 5000 atoms. The statistical 
accuracy in the lifetime determination reaches 20%. Data collected in 2002 allow us to increase 
the statistical accuracy up to 15%. 
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Prospects for the study of the Kf A decays at the NA48/2 

experiment 

C. Cheshkov and G. Marel 

DSM/DAPNIA - CEA Saclay, F-91191 Gif-sur-Yvette, France 

It is well known that the study of the K ± — > ■K + n~e ± i' e decay parameters can provide an 
important information about the 7nr interaction near threshold. Besides the extraction of the 
decay form-factors which are valuable input for xPT, the measurement of the 7i7i scattering 
length dp is recognized as a crucial cross-check with the current understanding of chiral sym- 
metry breaking of QCD and as an accurate estimate of the size of quark condensate < 0|mt|0 > 
[1]. The two most significant experiments so far were carried out by the Geneva-Saclay [2] and 
E865 collaborations [3] and are based on 3 x 10 4 and 4 x 10 5 collected — > 7r + 7r _ e ± f e events. 
Since the present experimental uncertainty on a° is significantly larger than the theoretical one 
it is of great interest to acquire an additional high statistics and good quality experimental 
data in this decay mode. 

The NA48/2 experiment is approved as an extension of the NA48 experiment at the CERN 
SPS. Its main goals are the search for direct CP violation in the Dalitz plot asymmetries in 
the decays of opposite charged kaons, the measurement of the scattering length a° and the 
analysis of various rare kaon decays. About 1 x 10 11 K + and 0.6 x 10 11 K~ decays are expected 
for around 100 days of data taking in 2003. The NA48/2 experiment is based on a slightly 
modified NA48 set-up. The kaon decays volume is contained in a 114m long vacuum tank 
which is followed by the NA48 detectors. The charged particles are detected and reconstructed 
by a high resolution magnetic spectrometer. It consists of four drift chambers and a dipole 
magnet with horizontal momentum kick of 120MeV/c. The corresponding track momentum 
resolution is given by ap/P[%] 0.48 ©0.02 x P[GeV/c\. The energy, position and time of the 
showers produced by the charged particles passed through the spectrometer are measured by a 
quasi- homogeneous liquid krypton electromagnetic calorimeter (LKr). Its transverse structure 
of p^1300 readout cells each with size of 2 x 2cm 2 is formed by electrodes extended with an 
accordion geometry. LKr has a longitudinal size of 127cm (equivalent to 27X and ps 2A) 
and a projective tower geometry pointing to the decay volume. Its energy resolution for e/m 
showers is cr E /E = 0.42%©3.2%/ v / £©0.09/£ (E in GeV). A detailed description of the whole 
NA48 set-up can be found elsewhere [4]. 

In addition to the present NA48 apparatus, two essential new elements will be introduced to 
the NA48/2 set-up [5]. The first one is an achromatic beam-line capable to transport and focus 
simultaneously 60GeV/c±10% K + and K~ beams. A kaon beam spectrometer (KABES) will 
measure the kaon beam sign, momentum and direction in order to improve the reconstruction 
of the Kf A and to recover the — > 3n with one pion escaping detection [5]. KABES consists 
of three so-called stations - two in the beam-line achromat and one downstream. Each station 
is formed by two time projection MICROMEGAS chambers [6] with opposite drift directions. 
The detector prototype has been successfully tested in July 2002 in the environment of an 
extremely high particle rate of 20-30MHz (comparable to that of the future kaon beam-line) 
showing excellent time (<lns) and space (<80/im) resolutions. 

Because of the small branching ratio of K ± — > n + 'K~e ± v e an important point in the 
systematic free measurement of a° is the reduction of the major background coming from 
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Table 1: Resolution on Cabibbo-Maksymowicz 
variables. 



K ± - 
K ± - 
K± -> 
K±- 
K ± - 
K ± - 



— ► 7T 7T 7T~ 
• 7T + 7r _ 7T ± - 
> T^TT - 
TT^TT ?! — i 

e i/7T° — > 



(r) 

e + e~7 
-> e + e~7 
> e + e~7 
> e + e~7 



Background 



< 1% 

< 0.4% 

< 0.3% 

< 0.1% 

< 0.1% 

< 0.1% 

< 0.2% 



Table 2: Background contaminations. 



— > 7r + 7r 7T 111 decays with 7T ± misidentified as an electron. The situation is even more com- 
plicated since the background populates the low ttti invariant mass region which is the most 
sensitive to 0°. To cope with the relatively large expected — > tv + 7v~tv ± background the 
NA48/2 collaboration has introduced a Neural Network (NN) based approach [7] of exploit- 
ing the whole available LKr information about the longitudinal and transverse electromagnetic 
and hadronic showers development. Each electron candidate is described by the charged track 
measured by the spectrometer and the associated LKr shower parameters and is passed then 
as an input to a specially designed multilayer NN. The candidate is identified as an electron or 
rejected as a pion taking into account the probability of electron hypothesis given by NN out- 
put. The NN was developed, trained and tested using pure experimental samples of electrons 
and charged pions. By applying the described above method a pion rejection factor of 3500 at 
94% electron detection efficiency is achieved. 

In order to estimate the expected acceptance, resolution and background for Kf A decays a 
GEANT based Monte-Carlo simulation of the apparatus has been used. The K ± — > ir + 7r~e ± i> e 
decay has been described by three real form-factors (g = g/ f s , g' — g'/f s , h = h/f s ) and one 
phase (5 = Sq — S\). The values for these parameters were taken from [2]. The scattering 
length a[j has been extracted following the procedure in [2]. The selection of K — > n+'K^e^Ve 
events has been made by applying a set of standard acceptance and kinematic cuts (including 
a cut on the transverse momentum of the three charged tracks) as well as the NN electron 
identification. The expected acceptance of ~ 30% has been found relatively flat upon the 
Cabibbo-Maksymowicz variables leading to a conclusion of small systematic effects due to the 
inexact knowledge of the detector response. Taking into account the expected charged kaon 
flux in 2003 the expected K ± — > n + iT~e ± b , e decay statistics is > 1 x 10 6 and the corresponding 
statistical uncertainty on a° is less than 0.01. The estimated resolution and background con- 
taminations are presented in Table ^ and Table respectively. As one can see the extremely 
good performance of the charged particles spectrometer reflects in the high resolution on in- 
variant mass variables. On the other hand, the relatively high beam momentum of 60GeV/c 
leads to moderate resolutions upon the angle variables which are nevertheless comparable to 
those in the previous low energy experiments. Due to the high e/ir separation provided by the 
NN method the — ► 7i + ii^ii ± background is expected to be less than 1%. A detailed analysis 
of this background has shown that the corresponding systematic effect in the measurement of 
an is less then 0.01. 
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Figure 7: The 3n invariant mass distributions for the data (solid line) and MC (stars). 



In the framework of the preparation for the future experiment, in September 2001 the 
NA48/2 collaboration has performed a 4 hours charged kaons beam test run. The already 
tuned K eA decays selection and reconstruction has been applied to the collected during this 
run data. 1477 events passed all the selection criteria including the NN electron identification 
procedure. The corresponding branching ratio has been evaluated using — > 7r + 7r 7T decay 
as a normalization channel. The obtained preliminary result of BR= (3.75 ± 0.25) x 10 -5 is in 
good agreement with the world average value [8]. A careful study of the differences between 
the data and Monte-Carlo M 37r distributions has confirmed that the — > 3tt background is 
less than 1% (Figured). 

To summarize, the NA48/2 experiment is planning to collect > 1 x 10 6 K ± — > 7i + 7i^e ± i > e 
decays for one year data taking in 2003 and to measure the irn scattering length a[j with a statis- 
tical precision of <0.01. The excellent e/ir separation power combined with the good detector 
resolution will allow to keep and to control the systematic uncertainty to the level of less than 
0.01. Thus the future measurement of a[] by the NA48/2 together with the forthcoming results 
from the DIRAC experiment [9] are expected to improve significantly the present experimental 
picture on this important parameter. 
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Future experimental investigation of the tt + 7T atom 

L. Nemenov 

CERN, CH-1211, Geneva 23, Switzerland and JINR, Dubna, Russia 

The difference AE n between the energy of atomic ns and np states of A 27T includes the 
contributions of the vacuum polarisation AE™ , and of the strong interaction AE^ effects, 
where AE^ ~ 2ao + a 2 , and ao and a 2 are the s-wave nn scattering lengths with isotopic-spin 
quantum numbers 0,2 [1,2]. The value of AE'^ ac is well known from QED calculations [1,3]. 
For this reason the measurement of AE n will give, in a model independent way, the value of 
2a + a 2 [4,5]. 

Some possibilities are discussed for measuring the np — ns energy splitting in (tt + tt~) atoms 
in experiments with relativistic particle beams, using the electromagnetic field control of an- 
nihilation process [6,7]. Significant resonant enhancement of annihilation probability may be 
observed in an oscillating field. The resonance provides much more spectacular changes in the 
decay rate than corresponding decay rate changes in a steady field. The position of resonance 
on the frequency axis gives an additional and more accurate information on the np-ns transition 
energy, than the data on the ratio between the field-free and dc-field-induced annihilation rates. 
Numerical results are presented for the states n = 2 and n = 3. 

[1] A. Gashi et al., Nucl. Phys. A 628 (1999) 101. 

[2] G. Efimov et al., Sov. J. Nucl. Phys. 44 (1986) 296. 

[3] A. Karimkhodzhaev and R.N. Faustov, Sov. J. Nucl. Phys. 29 (1979) 232. 

[4] G. Colangelo, J. Gasser and H. Leutwyler, Phys. Lett. B 488 (2000) 261. 

[5] M. Knecht et al., Nucl. Phys. B 457 (1995) 513. 

[6] L.L. Nemenov and V.D. Ovsiannikov, Phys. Lett. B 514 (2001) 247. 

[7] L.L. Nemenov, V.D. Ovsiannikov, E.V. Tchaplyguine, Nucl. Phys. A 710 (2002) 303-328. 



26 



Relevance of Kn atom measurements 



H. Sazdjian 

Groupe de Physique Theorique, Institut de Physique Nucleaire, 
Universite Paris XI, F-91406 Orsay Cedex, France 



It is pointed out that measurements of the lifetime and of the 2s — 2p energy level splitting 
may provide information about the combinations (a^ 2 — a^ 2 ) and (2a^ 2 + a^ 2 )/3 of the strong 
interaction S-wave scattering lengths, with isospins 1/2 and 3/2, respectively. Corrections to 
the nonrelativistic formulas of these quantities can be estimated with the use of the chiral 
effective lagrangian in the presence of electromagnetism. The combination {a\j 2 — a^ 2 ) is 
not sensitive at leading order to the chiral symmetry breaking mechanism, but may provide 
information about the low energy constant L 5 . The combination (2a^ 2 + a^ 2 )/3 is sensitive to 
the quark condensate value and to Zweig rule violating effects through its dependences upon the 
low energy constants (2L 6 + L 8 ) and L 4 . The electromagnetic and isospin breaking corrections 
are estimated to be of the order of a few percent in the lifetime value. For the 2s — 2p energy 
level splitting, the first main correction comes from vacuum polarization with a 25% effect. 
The next corrections come from the other electromagnetic and isospin breaking effects and are 
of the order of a few percent. 
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The Roy equations for the ttK system 

P. Biittiker 1 , S. Descotes 2 , and B. Moussallam 1 

1 Institut de Physique Nucleaire, Universite de Paris-Sud, F-91406 Orsay 
2 Laboratoire de Physique Theorique Hautes Energies, Universite de Paris-Sud, 

F- 91 406 Orsay 

In the past, nK scattering attracted the interest of experimentalists as well as theorists. Un- 
fortunately, the theoretical analyses of this process often suffered from the low statistics of the 
experiments and sometimes from over-simplified theoretical assumptions. On the other hand 
high statistics experimental data became available in the late seventies and eighties which have 
never been subject of a stringent theoretical analysis. Recently, there has been a revival in 
the interest in nK scattering: there are indications for a flavour dependence of the size of the 
quark condensate [1], and as the 7rK system is the most simple 5'?7(3)-process involving kaons, 
this process is a suitable place to test this dependence; effects of isospin violation in the above 
process are studied in the framework of ChPT [2] in order to obtain reliable predictions for 
the two S'-wave scattering lengths from the nK bound state experiment of DIRAC [3] ; the ttK 
fT-term has been analyzed [4] which can be regarded as an intermediate step towards ambitious 
calculations for the KN system; and finally nK scattering data were used to investigate the 
existence of a k resonance [5]. 

However, all these calculations strongly depend on reliable experimental data at very low 
energies. Unfortunately, the available high statistics data start only at about 800 MeV so that 
one is left with a gap at low energies. It has been shown for the nn system that solving the 
Roy-equations is a suitable tool to close this gap [6]. Here we apply this method to the ttK 
system. 

Assuming isospin to be a conserved quantity, nK scattering can be described in terms of the 
isospin-even and -odd amplitudes F + (s, t) and F~(s, t), respectively, the crossed channel, nn — > 
KK is given by F It=0 (s,t) and F It=1 (s,t). By combining fixed-t and hyperbolic dispersion 
relations with at most two subtractions for each of these amplitudes, one can derive twice- 
subtracted dispersion relations where the two S*-wave scattering lengths a$ and a$ play the 
role of the unknown subtraction constants. Projecting these dispersion relations onto partial 
waves, together with unitarity, yields a set of coupled integral equations for the s- and t-channel 
partial waves (Roy-equations) with two a priori free parameters a£ and [7], e.g. 

2V2m + a t r°° dt'\m.g\(t') 
3(m 2 + -m 2 _) + n km* T' D - t 

+ - / ds>{G w (t,s')lmf (s') + G u (t,s')Imf-(s')} + d{(t), 

where the G~- are the kernel functions and d\ is a so-called driving term accounting for the 
contributions of the higher partial waves. 

Using the available experimental data for the s- and t-channel above 1 GeV and the driving 
terms as input, one can show — in analogy to the 7T7r-system, see [6] and references therein - 
that the system of Roy-equations has a unique solution below 1 GeV. 

The equations then can be solved in several steps: using a Schenk-type parametrization 
for the s-channel waves and due to extended unitarity the t-channel partial waves can be 



9l(t) = 
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solved by applying Omnes-Muskhelishvili methods. These solutions for the t-channel waves 
then are used in s-channel equations and by tuning the parameters of the s-channel waves one 
tries to find the best input /output agreement for the Roy-equations. This procedure generates 
approximative solutions of the nK Roy-equations for some ranges of values for Oq and Oq ■ 
However, most of the solutions are not physical and preliminary results indicate that in 7r.fr 
scattering the Roy-equations do not support a universal band, which is a relation between the 
two ,S-wave scattering lengths, but one definite point in the a^-a^ -plane. Furthermore, these 
results indicate a value for ag = {a]J 2 — aj^ 2 )/3 which is somewhat larger than what has been 
predicted by S77(3)-ChPT [8], calling thereby for an independent experimental check of the 
value of Gig by the irK bound state experiment of the DIRAC collaboration [3] . 
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At low energies, QCD is dominated by the spontaneous breakdown of chiral symmetry which 
leads to the appearance of pseudoscalar Goldstone bosons ir, K, r\. Due to the mass hierarchy 
m 'u,d m s <C Ah, we can consider two different chiral limits: 

Nf = 2 : m u , irid — > , m s physical , Nf = 3 : m u , m^, m s — > . (9) 

In both limits, chiral order parameters can be defined, such as the quark condensate: H(Nf) = 
\im Nf (0|ww|0), and the pseudoscalar decay constant: F(Nf) = lim^.F^. 

The patterns of chiral symmetry breaking reflected by these parameters can be quite different 
in the Nf = 2 and Nf = 3 limits [1]. For instance, the condensates are related through: 
S(2) = S(3) +m s Z sca i ar , where Z sca \ aT > describes the violation of the Zweig rule in the scalar 
sector ++ . There have been recently hints of a significant decrease from S(2) to S(3), which 
would lead to a different low-energy behaviour of the Nf = 2 {jnx scattering) and Nf = 3 (7r.fr 
scattering) sectors. 

Because of the possibility of a small three-flavour condensate S(3), more care is needed 
when dealing with chiral expansions obtained from Nf = 3 chiral perturbation theory, e.g.: 

F 2 M 2 = 2mS(3)+m;K Q log(M Q / /U ) + c«( /U )] + F 2 MX, (10) 
F 2 M| = (m + m s )S(3)+m 2 ^ Q log(M Q / yU ) + 4(/i)] + F 2 M|^. (11) 

Next-to- leading-order (NLO) terms contain chiral logarithms (Q = tt, K, 77), and two low-energy 
constants L 8 and L 6 (in cp). L e is related to Z scallir and measures how badly the Zweig rule is 
violated in the scalar sector. The remainders dp collect NNLO and higher-order terms. 

We want to take into account the possibility of a suppressed three-flavour condensate S(3). 
We will assume an overall convergence of the chiral series [NNLO remainders d n , dx ~ (30%) 2 ~ 
10%], but we do not make assumptions about the relative sizes of leading- and next-to-leading- 
order terms. We simply consider the chiral series (fTU|) and (jllj) as exact identities. 

We can first express the three-flavour quark condensate in terms of F(3), r and L 6 [1]: 

2mS(3) 1 - e(r) - d r - M M " (r , 2) 1 = fC?j = * it 

where r = m s /m, Z(3) = [F(3)/F 7T ] 2 , e(r) < 0.2 for r > 15 and Lg rit depends (weakly) on r only 
(Lf\M p ) = -0.26 • 10~ 3 for r = 25). The NNLO remainder d is estimated as d ~ d„ ~ 10%. 
The square root arising in eq. (|T2*|) indicates the nonperturbative resummation of Zweig-rule 
violating effects. It is expanded in the standard treatment of xPT, assuming £ 1 - inspired 
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by large- N c expectations. But a small shift of L 6 from Lg nt of a few 10~ 3 would be magnified 
by the huge scaling factor 64M^(r + 2)/F% ~ 4 ■ 10 3 , leading to a damping of X(3). 

For the two-flavour quark condensate, the chiral expansions of F%M% and F^M\ yield [2]: 

X(2)[l - 4] = X(3) + -^{[1 - X(3)] - e(r) - d - /[X(3)/Z(3)] 2 } , (13) 

where d^ ~ ± 3% and / is a function of chiral logarithms in the Nf = 2 chiral limit 
(0.03 < / < 0.06). The two-flavour quark condensate X(2) is therefore the sum of two terms: 
a genuine condensate X(3) and an induced condensate proportional to m s . For large r (r > 15), 
X(2) remains close to 1 for any size of X(3) due to a compensation from the induced condensate. 
This phenomenon can be interpreted in terms of the density and fluctuation of small eigenvalues 
of the Euclidean Dirac operator [2]. 

X(2) is weakly correlated with the pattern of Nf = 3 chiral symmetry breaking, but strongly 
correlated with the value of the quark mass ratio r [1], as seen from the plot on the left. 




r=ms/m r=ms/m 



Accurate measurements of irn scattering would pin down Nf = 2 chiral order parameters, 
but would only constrain mildly the pattern of Nf = 3 chiral symmetry breaking. The analysis 
of available tttt data [3] leads to X(2) = 0.81 ± 0.07, implying that r > 15. The Bayesian 
approach allows one to express more quantitatively this statement and to combine results 
concerning different observables [4]. The plot on the right side is a typical outcome, exploiting 
the 7T7T results for X(2) and Z(2) to determine the posterior probability of the quark mass ratio 
r. We see that the current irn data rule out small values of r (< 15) and favours slightly r ~ 20. 
Including observables from the Nf = 3 sector (ttK scattering, r\ — > 37r) within this framework 
should finally elucidate the Nf = 2 and Nf = 3 chiral structures of QCD vacuum. 
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interaction of 7r + 7r~-atom with atoms of matter 5 

A.Tarasov and O.Voskresenskaya 
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The Glauber theory for interaction of hydrogen-like elementary atoms (EA) with target 
atoms (TA), developed in the papers [1-3], is essentially based on the assumption that the 
Coulomb potential, created by TA, does not change during the EA-TA interaction. In other 
words, in this approximation all possible excitations of TA in intermediate and/or in finale 
states are completely neglected. 

The Glauber theory corrected to account these effects is too cumbersome. Here we list the 
simplest results of this theory concerning the total cross sections of EA-TA interactions. 



o M (i) = o&(i) + o% !0h (i) (14) 

vlo h ,incoh(i) = J ^r\%(r)\ 2 <fbF{b, s) coh , incoh (15) 

r(6, Z) coh = 1 - 2 cos [A X (b, s)} exp [-$(6, s)/2\ + exp [-$(6, s)\ (16) 

r(6, s)incoh = 1 - exp [-$(6, s)\ (17) 

Ax(6, 8) = ^J^ (e^+ _ e «K_) [i _ Sl (q)} (18) 

Hb, s) = ^rj 5^ (e*& - e^~) (e~^ + - e*) W(q u q 2 ) (19) 

Wfafa) = - q 2 ) - S^S^fc) + (Z - l)[S 2 (q u q 2 ) - S^S^fc)] (20) 

Si(q) = J d 3 re^ Pl (r) , J d 3 r Pl (f) = 1 (21) 

S 2 (qi, q 2 ) = J cPr^ne^^pfa, q 2 ) , f dh^r^) = p 1 (f 1 ) (22) 

b± = b±s/2, s = f ± (23) 



W(q,q)=S mcoh {q) (24) 

Here pi )2 are one-particle and two-particle electron densities of the target atom, cr^ hincoh are 
the total cross sections of EA-TA interactions without or with excitation of the target atom. 
$ in the above equations accounts the target atom evitation both in intermediate and in final 

5 The report was not presented because of a visa problem 
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states. If one put $ = then Eqs.(l-4) turns to the corresponding relations of papers [1-3]. In 
particular in this limit (y incoh = 0. Relative corrections to cr*°^ caused by including intermediate 
incoherent effects are of order 



and can be successfully neglected. The same is also true for all partial coherent cross sec- 
tions. From these estimations it follows that the theory of papers [1-3] provides quite accurate 
description of the coherent sector of EA-TA interactions. 

As for incoherent interactions, it follows from Eq.(4), that they can be described by the 
Born approximation with the relative accuracy of order 



The results of performed analysis can be summarized as follows: 

1. for description of the coherent interactions of EA with TA it is enough to use the simplified 
version of the Glauber theory neglecting effect of intermediate excitation of TA; 

2. for description of the in coherent interactions of EA with TA it is enough to use the Born 
approximation. Just such prescriptions based on an intuitive consideration have been 
proposed by authors of paper [3]. 
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Concluding remarks 



H. Leutwyler 

Institute for Theoretical Physics, University of Bern, Sidlerstr. 5, CH-3012 Bern, Switzerland 

In the first part of the talk, I illustrated the relevance of low energy pion physics for our 
current understanding of the basic interactions, using the example of the Standard Model 
prediction for the muon magnetic moment. I then discussed the remarkable recent progress 
achieved in mi scattering, which now has become a precision science: in combination with 
the low energy theorems of Chiral Perturbation Theory for the S-wave scattering lengths, the 
dispersion relations for the partial waves (Roy equations) pin down the scattering amplitude 
to an amazing degree of accuracy [1-3] - for once in strong interaction physics, theory is ahead 
of experiment. The recent K ei data from Brookhaven [4] offer a rough test of one of the 
hypotheses that underly the theoretical framework, as they show that the quark condensate 
indeed represents the leading order parameter of the spontaneously broken chiral symmetry. 
Generalized ChPT has served its purpose and can now be dismissed. 

The Roy equations involve two subtraction constants, which may be identified with the two 
S-wave scattering lengths. K e4 decays are sensitive to one particular combination of these. The 
forthcoming results from NA48/2 are eagerly awaited as they will narrow down the uncertainty 
in the experimental information about this combination [5]. The low energy theorem for the 
scalar radius involves a different combination, but this theorem has not been subject to exper- 
imental test. The importance of such a test was thoroughly discussed at this meeting [6-8]. A 
somewhat more precise determination of the I = 2 7T7t phase shift may become possible at the 
COMPASS detector at CERN [9] and, as pointed out by L. Montanet, there are beautiful data 
on other processes with two-pion configurations in the final state, such as D — > 3it, which may 
also help to narrow down the experimental uncertainties. 

The ideal laboratory for exploring the low energy properties of the pions, however, is the 
atom consisting of a pair of charged pions, also referred to as pionium. The DIRAC collabora- 
tion at CERN has demonstrated that it is possible to generate such atoms and to measure their 
lifetime [10]. Since the physics of the bound state is well understood [11-13], the low energy 
properties of the scattering amplitude can unambiguously be determined by studying these. It 
would be most deplorable if this beautiful project were aborted before its physics potential [14] 
is tapped. In particular, pionium level splittings would offer a clean and direct measurement 
of the second subtraction constant. Data on nK atoms would also be very valuable, as they 
would allow us to explore the role played by the strange quarks in the QCD vacuum [15,16]. 

Experimentally, pionic hydrogen has been explored very successfully [17]. Also, there is 
significant progress in the physics of the bound state [18-20]. The main problem here is that 
the ChPT predictions for the scattering lengths have large uncertainties. At this time, accurate 
results can only be obtained in the unphysical region of irN scattering. Dispersion theory is 
needed to confront these with data at or above threshold. An analogue of the Roy equations 
has been formulated [21], but the framework yet needs to be implemented. A comprehensive 
dispersive analysis of the irN scattering amplitude was carried out by G. Holder and collab- 
orators 20 years ago [22]. An update of that analysis is under way, but will take some time 
for completion [23]. Currently, the Standard Model predictions for pionic hydrogen are too 
uncertain to draw physics conclusions from the available, beautiful data, concerning the size of 
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the (T-term, for instance. The same applies to the predictions for KN atoms, which are under 
investigation in the DEAR experiment at Frascati [24]. 

There is plenty of work to be done in the field of hadronic atoms. Quite a few experimen- 
tally accessible quantities can be calculated from first principles. By confronting theory with 
experiment, we can critically examine our understanding of the laws of nature. 
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